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We compute cosmic microwave background angular power 
spectra for scaling seed models of structure formation. A 
generic parameterization of the energy momentum tensor of 
the seeds is employed. We concentrate on two regions of pa¬ 
rameter space inspired by global topological defects: 0(4) 
texture models and the large-A limit of 0(N) models. We 
use \ 2 fitting to compare these models to recent flat-band 
power measurements of the cosmic microwave background. 
Only scalar perturbations are considered. 

PACS: 98.80-k, 98.80Hw, 98.80C 

Inflation and topological defects are two families of 
models to explain the origin of large scale structure in 
the universe. In models with topological defects or other 
types of seeds, fluctuations are generated continuously 
and evolve according to inhomogeneous linear perturba¬ 
tion equations. Seeds are any non-uniformly distributed 
form of energy, which contributes only a small fraction 
to the total energy density of the universe and which in¬ 
teracts with the cosmic fluid only gravitationally. We 
are particularly interested in global topological defects 
playing the role of seeds. 

Cosmic microwave background (CMB) anisotropies 
provide a link between theoretical predictions and obser¬ 
vational data, which may allow us to distinguish between 
inflationary perturbations and defects, by purely linear 
analysis. On large angular scales, both families of mod¬ 
els predict an approximately scale-invariant Harrison- 
Zel’dovich spectrum [|i]j^]. For inflationary models this 
can be seen analytically. Scale-invariance for defects was 
discovered numerically 11 ; simple analytical arguments 
are given in [||. At small angular scales (0°2 ~ 6 ~ 1°), 
the predictions of inflation and defects are different. 
CMB observations at these scales may soon be sensitive 
enough to distinguish the two families of models. 

In a recent work [Q , two of us investigated the general 
behavior of CMB anisotropies induced by seeds. Here, 
for simplicity, we restrict ourselves to scalar type pertur¬ 
bations. Thus, the models presented in this work are not 
close approximations to the 0(4) texture model for which 
the Sachs-Wolfe (SW) plateau is dominated by vector and 
tensor modes jsflij. 

This restriction does not render our work uninterest¬ 
ing. There may very well be models with scaling seeds 
leading to small vector and tensor contributions, e.g., 


due to symmetry constraints (spherical symmetry) or in 
models with non-relativistic seeds. Here, we assume a 
completely phenomenological standpoint: we investigate 
whether models with purely scalar seeds can reproduce 
the data. In subsequent work we plan to study how 
severely vector and tensor contributions are restricted by 
the data. However, the models already excluded on the 
basis of our work, will not be resurrected once vector and 
tensor modes are included. 

In our models, we characterize the energy momentum 
tensor of the source by four seed functions which we term 
f p , f p , f v , and A, defined by (see @0) 


= M 2 f p 

(1) 

= M 2 f Vii 

(2) 

= M 2 [{f p — (l/3)Af n }'yij + fn\ij\ , 

(3) 


where A denotes the Laplacian and | is the covariant 
derivative with respect to the metric 7 of three space. 
M is a typical “mass”, or energy scale, of the seeds. The 
gravitational strength of the seeds is characterized by the 
dimensionless parameter e = 4ttGM 2 . The superscript 
indicates that only the scalar contribution to 0,;o and 
Qij is included here. Since seeds interact with other mat¬ 
ter components only gravitationally, the seed functions 
satisfy the covariant conservation equations jn| 

fp - A/„ + (a/a)(f p + 3 f p ) = 0 (4) 

fv + 2(d/o)/„ -f p - (2/3)A/ 7r = 0 , (5) 

where a is the scale factor and dot stands for derivative 
with respect to conformal time t. 

We define “scaling seeds” to be seeds for which the 
power spectra (|/.| 2 )(fc, t) are, up to an overall power of 
t determined by dimensional reasons, functions of x = kt 
only. Thus, the power spectra of the functions /. are of 
the form 

(\f P \ 2 ) = t ~ 1 F?(x) 

(\f P \ 2 )=t- 1 F 2 (x) 

(| f v \ 2 )=tFi(x) 

(\U\ 2 )=t 3 F 2 ( x ) . (6) 

Furthermore, we require that the seeds decay on sub¬ 
horizon scales. This behavior is found in simulations for 
the seed functions of global textures and is also supported 
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by the large- N limit of global O(N) models [p~2[[l~3[ . Sim¬ 
ple analytical arguments indicate that all types of models 
with scaling seed functions which decay fast enough in¬ 
side the horizon lead to a Harrison-Zel’dovich spectrum 



Numerical simulations of global O(N) models show 
that on super-horizon scales (x <C 1 ), Qij and ©oo have 
white noise spectra , whereas |©io | 2 behaves like k 2 . Fur¬ 
thermore, the power spectra of the functions /. do not 
depend on the direction of k. Thus the F,'s are even 
functions of x = kt. Consequently, iq, iq, F 3 —► const., 
while Fi oc 1/x 2 for x —» 0. Since the energy momentum 
tensor of the seeds decays inside the horizon, Fi —> 0 
for x —y 00 . In this work we approximate the random 
variables /. by the square root of their power spectra. 
Motivated by numerical simulations and the considera¬ 
tions described above we model the functions F\ , F 3 as 


F! 


A 1 

1 + aix ni 


F 3 


a 3 

1 + a 3 x ns 


(7) 


F -2 and F 4 are then given by energy momentum conser¬ 
vation, Eqs. ([|, ||). 

The gravitational action of the seeds is determined by 
the induced Bardeen potentials, which are not only due 
to the seeds but also contain contributions from the mat¬ 
ter and radiation fluids. Once the fluid perturbations 
and the Bardeen potentials are determined, one calcu¬ 
lates CMB anisotropies by standard methods for each 
model parametrized by (Ai, aq, n\, A 3 , a 3 , n 3 ). For de¬ 
tails see Ref. |7]]. 

In this letter we present the results of a parameter 
study and we fit to the observational data available. We 
compare the anisotropy power spectra obtained in our 
models with observations. The cosmological parameters 
used are h = 0.5, fls = 0.0125 h ~ 2 ,12 = 1 and A = 0. We 
are thus considering scaling seed models in the context 
of flat cold dark matter universes. 

We investigate two types of models. In the first one we 
choose n\ = 2, 713 = 4; a choice supported by numerical 
simulations of global textures [Q; we refer to this first 
type of scaling seed models as model A. In the second 
one we set n 1 = 5/2 and n 3 = 7/2; a result obtained 
analytically in the large- N limit of O(N) models [ ^2| . 
We call it model B. 

We set the arbitrary normalization by fixing A\ = 1 
and we vary A 3 i= A 3 /A\). To make the calculations 
feasible, we further reduce the remaining 3-dimensional 
parameter space (A 3 , aq, a 3 ) to 2-dimensions (A 3 , aq), by 
setting a 3 = oq/2. Thus, our fits are displayed as contour 
plots in the oq — A 3 plane (see Figs. 3 and 4). We have 
also investigated a 3 = 2 aq and obtained qualitatively 
similar results. 

The A 3 dependence is similar for both types of models. 
For A 3 ~ 0.06, the relative amplitude of the acoustic 


peak, with respect to the SW plateau, decreases as A 3 
decreases (see Figs, la and c). 

There is a particular value of the constant A 3 for which 
the coefficient of the 1/x 2 term in F 4 vanishes. For x <C 
1, one obtains from Eqs. ([|,|]) F 4 ss A^/x 2 with A 4 = 
— (3/8)Ai(l + I 8 A 3 /A 1 ) in the matter dominated era. If 
A 3 = (—l/18)Ai = A crit ~ —0.06, we thus find A 4 = 0 
[|7j. This 1/x 2 term dominates on super-horizon scales 
and (for big enough aq) determines the amplitude of the 
SW plateau. Its absence is thus expected to lead to a 
higher relative amplitude of the first acoustic peak, which 
is well visible in Figs, la and lc. 



FIG. 1. Parameter dependence of the calculated power 
spectra. To illustrate the dependence on A 3 , we choose val¬ 
ues around A cr it = —0.06 (see text). The dotted line is a 
polynomial fit to the data and is the same in each panel as a 
reference. 

In our exploration of parameter space we vary —1.0 < 
A 3 < +1.0 for model A and —1.2 < A 3 < +0.5 for model 
B. For oq we choose the parameter range 0.01 < oq < 2.0 
for model A and 0.001 < oq < 0.141 for model B. We 
normalize the power spectra at £ = 10 by fitting to the 
data. 

We use x 2 fitting to compare our models to recent 
flat-band power measurements of the CMB. The method 
and a compilation of the data is described in detail in 
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FIG. 2. Flat-band power measurements used in this anal¬ 
ysis. The best-fit models of two scaling seed models are also 
shown: one most relevant to model A (Fig. 3) and the other 
for model B (Fig. 4). 

p|| . The most recent data and improvements to the x 2 
method are described in jnjjnj]. 

Fig. 2 plots the data used along with the best-fit power 
spectra for the two types of models described above. 
These best-fit models are indicated by the “X” in Figs. 3 
and 4. Fig. 3 is a contour plot in the ax — A 3 plane for 
model A while Fig. 4 is the analogous one for model B. 

There are 32 data points and 28 degrees of freedom (28 
= 32-2 (fitted defect parameters) - 1 (normalization) - 1 
(sliding Saskatoon absolute calibration)). Model A yields 
a x 2 minimum value y 2 lin = 21.2, while model B yields 
Xmin = 27.1. Thus the fits are good. 

To interpret our y 2 fits, we first note the following: if 
perturbations decay fast enough, the height of the first 
acoustic peak is determined by / p +3/ p ~ Ax/A, while the 
SW plateau is fixed by A 4 . The \ 2 is then not expected 
to be very sensitive to a. This indeed is seen in model B 
for ol\ ~ 0.03 (Fig. 4 and also Fig. 1 b) and in model A 
for ax ~ 0.4 (Fig. 3.). 

We also have investigated the dark matter power spec¬ 
tra for some of the parameter space of our models. For 
the best fitting models we obtain tr 8 ~ 0.6 ± 0.2 which 
is in reasonable agreement with observations. Values of 
(A3, ax) which are excluded from the Ce analysis, often 
lead to far too small values of trg. We found, however, 
that the bend in the power spectrum of our “best mod¬ 
els” lies at somewhat smaller scales than in the analysis 
of APM and IRAS galaxies done by Peacock Q . Since 
his analysis assumes Gaussian statistics, we are very re¬ 
luctant to draw any strong conclusions from this compar¬ 


FIG. 3. x 2 fit of model A to the CMB data. The “X” marks 
the Xmin ; the grey area is the (Xmin + 1), ~ la-likelihood 
region. Contours marked “68”, “95” and “99” refer to good- 
ness-of-fit contours. For example, under the assumption that 
the errors on the flat-band power measurements are Gaus¬ 
sian, the probability of obtaining a y 2 value less than the 
value obtained on the “95” contour is 95%. 

ison. 

In the matter dominated era, numerical simulations 
of 0(4)-textures give Ax ~ 4,cti ~ 0.012 and A 3 ~ 
0.37, 0:3 ~ 0.05 |l^|. For this model, the scalar contri¬ 
bution to the SW plateau is ~ l.le 2 at t ~ 10 and the 
height of the first acoustic peak is about ~ 5e 2 . On the 
other hand, full simulations, which include vector and 
tensor modes jf|, lead to an amplitude of the SW con¬ 
tribution on the order of i(t + 1 )Ce ~ S^e 2 - Therefore 
the vector and tensor parts contribute more than ~ 80% 
to the SW plateau, while they are not expected to influ¬ 
ence the acoustic peaks. In the full texture models, the 
acoustic peaks are thus expected to be substantially too 
low to fit the data. This result was pioneered in || and 
has now been confirmed by full numerical simulations |5| . 
In Ref. H decoherence Jl^] has also been taken into ac¬ 
count, which further reduces the acoustic peaks without 
influencing substantially the SW plateau. 

The fact that we obtain a parameter range compatible 
with currently available observational data is clearly not 
in contradiction with the result of || and |J, namely 
that the acoustic peaks for the conventional 0(4)-texture 
model are very small or even completely absent. 

Fig. 3 shows that small negative values of A 3 (—0.1 < 
A 3 < 0) are preferred. For ax ~ 0.4 the result does 
not depend strongly on ax . For smaller values of ax, 
(ax ~ 0 . 2 ) somewhat larger values of |^ 31 are preferred 
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FIG. 4. Same as previous figure except here we are fitting 
model B to the data. 

since the acoustic peaks are already enhanced due to the 
slower decay of the seed functions. 

Fig. 4 is the model B analog of Fig. 3. Also here, a 
positive A 3 is excluded for aq ~ 0.01. A value of A3 ~ 
—0.1 is generically preferred. Note however, that the 
X 2 “landscape” within the parameter range explored in 
Figs. 3,4 is rather flat, and values within the “68” contour 
are reasonably compatible with current data. 

In this letter, our aim was not to test whether a given 
model with topological defects can fit the data. We 
wanted to investigate, whether present observations of 
CMB anisotropies can already rule out a generic class of 
seed perturbations constrained just by energy momen¬ 
tum conservation and scaling arguments. Our analysis 
indicates that the answer to this question is no. 

As a continuation of this work, we plan to include vec¬ 
tor and tensor perturbations as well as decoherence in 
our models. We also want to study whether there are 
more severe restrictions on defect models than just en¬ 
ergy momentum conservation and scaling; for example, 
to see whether the vector component always dominates 
the level of the SW plateau. 
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